# Figure S1: Alternative income threshold (unconditional AMCE)

fig3_inc2_dat = import(here("data","fig3_inc2_dat.rds")) |> 
  invert_group() |> fix_healthcare()

fig3_inc2 = plot_conjoint(conjoint = fig3_inc2_dat,
                          is_result = T,
                          lab = policies, 
                          weights = dta_grp$demog_weight[!is.na(dta_grp$inc2_grp)],
                          use="all",
                          title=NULL,
                          grouped_data=T) +
  labs(subtitle = "Income (>200K)", x = NULL, y = NULL)

print(fig3_inc2)

ggsave(here("figures","appendix","figure_s1.png"), fig3_inc2,
       dpi = 600, units = "in", width = 10, height = 9)

ccat("See 'figure_s1.png' for Figure S1 output")